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Abstract 



We consider the problem of choosing Euclidean points to maximize the sum of their weighted 
pairwise distances, when each point is constrained to a ball centered at the origin. We derive 
jy\'. a dual minimization problem and show strong duality holds (i.e., the resulting upper bound is 

tight) when some locally optimal configuration of points is afimely independent. We sketch a 
polynomial time algorithm for finding a near-optimal set of points. 



1 Introduction 

We consider the following maximization problem P(n,w,' 



(N ■ maximize^ >^ w(i,j)d(pi,pj) 

£5: ,. .. J Pi G IR"- 1 (* = l,..,n); 

^ . subject to < 

[ \\Pi\\ < i(%) (* = l,..,n). 

Here each w(i,j) > and each £(i) > is fixed, d(p,q) denotes the Euclidean distance between 
points p and q, and ||p|| denotes the Euclidean length (distance from the origin) of point p. 
H ' We derive the following dual problem D(n,w,£): 



minimize r x .i 
subject to 



y- w 2 (i,j) x 



E^ 2 ( 



l)Xi x 



E 



\ ±<i<3<n XiX 3 \ t=l \ i=l 

xi S M (i = 1, ..,n); 

xi > (i = 1, .., n). 
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Throughout the paper, jj is defined to be 0. 

We show that the value of the maximization problem is at most the value of the minimization 
problem. We use a physical interpretation of the two problems to show that the values are equal 
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provided the maximization problem admits a set of points {pi} that is both affinely independent 
and stationary (i.e., the gradient of the objective function is a nonnegative combination of the 
gradients of the active constraints, a necessary condition at any local maximizer of P(n,w,£)). 

We sketch how a near-optimal solution to the problem can be found in polynomial time via the 
ellipsoid method. 

2 Related Work 

The case w(i,j) = £(i) = 1 (in which the optimal points are given by the vertices of the regular 
n-simplex, achieving a value of nJ Q)) was previously considered by [3]. Our Lemma [T] generalizes 
a bound in that paper. 

Specific instances of P(n, w, £) were studied to obtain geometric inequalities that were used to 
analyze approximation algorithms for finding low-degree, low-weight spanning trees in Euclidean 
spaces [2]. 

Goemans and Williamson [lj consider related problems with applications to approximating the 
maximum cut in a graph and to maximizing the number of satisfied clauses in a CNF formula. We 
modify their approach to solving their problems to obtain a polynomial time algorithm for ours. 

3 A Dual Problem 

Lemma 1 For any n, w, and £, the value of the maximization problem P(n,w,£) is at most the 
value of the minimization problem D(n,w,£). 

Proof: Fix any n, w, and £. Fix any set of points {pt} and values {xj} meeting the constraints 
of P(n,w,£) and D(n,w,£), respectively. Let A(i,j) = ^j== and B(i,j) = ^/XiXjd(pi,pj) for 
1 < i < J ' < n - Then, by the Cauchy-Schwartz inequality A ■ B < \\A\\ x \\B\\ (where A and B are 
interpreted as (^-dimensional vectors, and • denotes the dot product): 



^w(i,j)d(pi,pj) < 

i<j 



\ 
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It remains only to show 



Expanding the left-hand side 



i<i \ % / \ % 



,2 (i)x,. 



^x i Xjd 2 (p i ,p j ) 



i<j 



- J2 XiXj{pi -Pj)- {pt - Pj) 

- Y^ XiXj(pi ■ pt - 2pi ■ pj + pj ■ pj) 



< ^2xiXj(£ 2 (i) -pi-pj) 

< (e^)x(e^ 2 w)- 



(2) 
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Lemma 2 Fix any n, w, and £. Suppose the maximization problem P(n,w,£) admits a set of 
points {pi} that is both stationary and affinely independent. Then the values of the two problems 
are equal. Further, there exists {x{\ such that 



%iPi 



E 



w{hJ. 



Pi - Pj 

d{pi,Pj) 



(where x% = in case \\pi\\ < £i, and w(i,j) = w(j,i) and w(i,i) 
optima for the two problems. 



(4) 



0), and {pi} and {xi} are global 



Proof: Fix any n, w, and £. Consider the objective function $({pj}) = J2ij w (^J)d(Pi,Pj) of 
P(n,w,£). That {pi} is stationary means that the gradient of <£ is a nonnegative combination of 
the gradients of the constraints of P(n,w,£) active at {pi}. By elementary calculus, the gradient 
of $ consists of a vector /j for each point pi, with each fi equal to the right-hand side of (JH). The 
only constraint on pi is \\pi\\ < £(i), whose gradient (again by elementary calculus) is a nonnegative 
multiple of p«. Thus, for each i, there exists an xi > such that ((4]) holds. Note that if ||pj|| < £(i), 
then the constraint is not active, so that fi must be the zero vector. In this case we take X{ = 0. 

We will show that each inequality in Lemma Q] is tight for these {pi} and {x^}. Inequality ([3]) 
is tight because, by (|3J), J2i x iPi is the zero vector. Inequality ([2]) is tight because \\pi\\ < £(i) only 
if Xi = 0. 

Inequality ([T|) is tight provided the vector A (in the proof of Lemma [1]) is a scalar multiple of 
B. Assume {p{\ is affinely independent. Then, considering {xi} and {pi} fixed and {w(i,j)} as the 
set of unknowns (i.e., reversing their roles), (Jl]) uniquely determines each w(i,j). Since 



w(i,j) 



XjXjd(pi,pj] 

J2k x k 



[1 < i < j < n) 



(5) 



is consistent with (jl|) (check this by substitution for w(i,j) in (JH), it follows that ([5]) necessarily 
holds. Thus, A is a scalar multiple of B and Inequality ([T]) is tight. D 

A physical model for the quantities involved is as follows. Consider a physical system of n 
points {pi}. Each point pi is constrained to a ball of radius £{i) centered at the origin. For each 
pair of points (pi,Pj), Pi repels pj (and vice versa) with a force of magnitude w(i,j). 

Under this interpretation, each vector fi in the proof corresponds to the force on pi, and Xi is 
the magnitude of this force, divided by \\pi\\. 



4 Solving P(n,w,£) in Polynomial Time 

If the instance of P(n, w, £) is small or has a high degree of symmetry, the dual problem D(n, w, £) 
might yield a function that can be minimized directly by symbolic methods. In general, it is 
possible to solve P(n, w, £) (to any given degree of precision) in polynomial time using semi-definite 
programming, following the approach in pQ. 

Those authors consider a related problem GW(w,n): 

maximize^} ^ w(i,j)d 2 (pi,pj) 

l<i<j<n 



subject to 



Pi £ R n (t = l,..,n); 

||pi|| = 1 (i = l,..,n). 



The authors show how to solve this problem in polynomial time by formulating it as a semi-definite 
program, and how to round a (near-)optimal set of points {pi} to obtain an approximate solution to 
a corresponding max-cut problem. This approach yielded the first polynomial-time approximation 
algorithm achieving a performance guarantee better than two for the max-cut problem. 

We briefly sketch their aproach for solving GW(w,n) and how it can be modified to solve 
P(w,n,£). The connection between sets of points and positive semi-definite matrices is the follow- 
ing: an n x n symmetric matrix Y is positive semi-definite if and only if there exists a set of n 
points {pi} in K n such that Yij = pi ■ pj. Thus, GW(w, n) is equivalent to following: 

maximizejy} y^w(i,j)(2 — 2Yjj) 

ij 

Y is an n x n symmetric, positive semi-definite matrix; 
subject to { 

Yn = 1 (i = l,..,n). 

The space of n x n symmetric, positive semi-definite matrices admits a polynomial time sep- 
aration oracle because a symmetric matrix Y is positive semi-definite if and only if x Yx > 
for each x G H n , and in fact it suffices to check each eigenvector x of Y. Thus, combining the 
constraint that Y is positive semi-definite with arbitrary linear inequalities on the elements of Y 
yields a convex space with a polynomial time separation oracle. Approximate feasibility of such a 
problem is testable in polynomial time via the ellipsoid method. Thus, GW(n) can be solved to 
near-optimality in polynomial time. 

A similar approach can be used to solve P(n,w,£) in polynomial time. In particular, P(n,w,£) 
corresponds to the following semi-definite program: 



maximizejy} J^ w(i,j)yYu + Yjj - 2Y^ 

Y is an n x n symmetric, positive semi-definite matrix; 
subject to < 

Yn < t(i) (i = l,..,n). 



Since Y^ij w(hJ)y/Yii + Yjj ~ 2Yij is a concave function in {Y^} whose gradient can be computed 
in polynomial time, the above program also admits a separation oracle sufficient to solve it to 
near-optimality in polynomial time using the ellipsoid method. 



5 Open Problems 

It would be interesting to obtain a more efficient algorithm for solving P(w, n, £) than is obtained 
by reducing to the ellipsoid method. Especially interesting would be a primal-dual algorithm along 
the lines of traditional "combinatorial" algorithms for solving or approximating linear programs. 
It is not clear how to achieve such algorithms in the semi-definite setting. 

Similarly, the only known method for achieving a better factor than two for the max-cut problem 
is by reduction to semi-definite programming. Goemans and Williamson leave open the problem of 
finding a more efficient algorithm that beats a factor of two. A more efficient algorithm for P(n, w, £) 
(with each £{i) = 1) would solve this, because applying their randomized rounding technique to 
P{n, w, €) also yields an approximation algorithm for max-cut with performance guarantee better 
than two. 

On the other hand, consider the generalization of GW(n, w) in which the objective function is 
replaced by Yl<ij w (h j)d 2+e (pi,Pj) for some e > 0. For e > 0, applying Goemans and Williamson's 
approach to this program rather than GW(n, w) would provide a better approximation to max-cut. 
Is the generalization solvable in polynomial time for some e > 0? 
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